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Foreword 

This  document  details  the  scientifie  progress  made  during  the  3-year  grant. 

Results  published  or  under  review 

1 .  Mathematical  modelling  shows  that  biogeochemical  feedbacks  in  microbial  communities 
ean  cause  oxie-anoxic  regime  shifts  (Bush  et  al,  in  collaboration  with  Profs  Huisman  and 
Muyzer,  University  of  Amsterdam,  manuscript  under  review  at  Nature  Communications) 

2.  Mathematical  modelling  shows  that  nonlinearities  in  microbial  population  dynamics  can 
cause  redox  regime  shifts  (Bush  et  al,  Biogeosciences  12,  3717-3724  (2015)) 

3.  Mieroeosm  experiments  show  that  initially  identical  replicate  communities  ean  diverge 
into  apparent  alternative  stable  states  (Pagaling  et  al,  manuseript  under  review  at 
Environmental  Mierobiology) 

4.  Mieroeosm  experiments  show  that  selection  history  affects  the  predictability  of 
community  development  (Pagaling  et  al,  ISME  J.  8,  19-30  (2014)). 

Results  still  to  be  published 

1 .  Experiments  show  the  time-dependent  development  of  mierobial  nutrient-cyeling 
eeosystems,  for  different  nutrient  conditions,  in  terms  of  both  community  structure  and 
chemical  composition. 

2.  Experiments  show  response  of  steady-state  mierobial  microcosm  communities  to 
environmental  change,  in  terms  of  eommunity  strueture  and  chemical  composition. 

3.  Experiments  show  response  of  microcosm  communities  to  addition  of  the  pollutant  PCP 
and  also  to  addition  of  the  PCP  biodegrading  organism  D.  hafniense  PCP-1. 

4.  Development  of  a  mathematical  model  to  mimic  the  effects  of  adding  organisms  that  can 
mediate  reductive  biodegradation  to  a  nutrient-eycling  mierobial  community. 

List  of  appendices,  illustrations  and  tables 

•  Eigure  1 .  Replicate  microcosms  created  from  freshwater  pond  sediment  and  water, 
supplemented  by  eellulose  and  sulfate. 

•  Eigure  2.  Analytical  predictions  of  a  ‘2  population’  model  for  a  biogeoehemical  cycle 
showing  that  it  predicts  redox  regime  shifts  as  a  funetion  of  the  reductive  and  oxidative 
driving  forces. 

•  Eigure  3.  Results  for  a  4-population,  2-box  mathematical  model  of  a  biogeoehemical 
cycle,  showing  that  this  model  also  predicts  redox  regime  shifts. 

•  Eigure  4;  Schematic  diagram  showing  the  mierobial  populations,  ehemical  species  and 
their  interactions,  in  our  mathematical  model  of  a  nutrient-cycling  ecosystem  in  which 
feedbacks  lead  to  alternate  stable  states. 

•  Eigure  5;  Dynamical  trajectories  generated  by  the  model  for  Eigure  4,  showing  that 
different  initial  system  states  can  produce  very  different  final  communities,  for  the  same 
set  of  parameters. 

•  Eigure  6;  Hysteresis  in  a  nutrient-cyeling  mierobial  eommunity.  (a):  Predietion  of  our 
mathematical  model  of  Figure  4;  as  oxygen  diffusivity  is  gradually  increased  the  system 
remains  for  a  long  time  in  an  anoxic  state;  as  oxygen  diffusivity  is  gradually  decreased 
starting  from  the  oxie  state,  the  system  remains  for  a  long  time  in  the  oxic  state,  (b);  Data 


gathered  from  Lake  Vechten  in  the  Netherlands,  showing  a  similar  hysteretic 
phenomenon. 

•  Figure  7.  Experimental  results  for  mieroeosm  strueture  as  a  function  of  time.  Summary 
plot  (a),  NMDS  plot  (b),  and  measure  of  community  diversity  (c),  showing  community 
composition  in  our  microcosm  communities  as  a  function  of  time,  as  determined  by  16S 
rRNA  gene  sequencing. 

•  Figure  8:  Experimental  results  for  steady-state  microcosms  with  varying  environmental 
parameters.  Diversity  (as  measured  by  Shannon  Index)  is  plotted  versus  relative 
abundance  in  our  set  of  microcosms  with  different  environmental  conditions,  for  the  three 
functional  groups:  photo  trophic  sulfur  bacteria,  sulfate  reducing  bacteria  and  oxygenic 
phototrophic  bacteria.  We  see  different  trends  for  the  different  groups,  suggesting  that 
their  community  structure  may  be  driven  by  different  ecological  mechanisms. 

•  Figure  9:  Schematic  illustration  of  our  mathematical  model  describing  a  nutrient  cycling 
microbial  ecosystem  supplemented  by  reductive  dehalogenators. 

•  Figure  10:  Experimental  results  for  microcosms  supplemented  with  PCP  and  D. 
hafniense:  the  persistence  of  D.  hafniense  in  the  microcosms  depends  on  the  level  of 
added  sulfate. 

•  Appendix  1 :  Equations  corresponding  to  our  4-population,  2-box  mathematical  model. 

•  Appendix  2:  Equations  used  to  describe  our  mathematical  model  for  the  nutrient-cycling 
microbial  ecosystem  with  alternate  stable  states. 

Statement  of  the  problem  studied 

The  bioremediation  of  contaminated  sites  relies  on  microbial  communities.  Different  populations 
of  microbes  undertake  different  chemical  transformations,  converting  potentially  harmful 
chemicals  via  a  series  of  intermediates,  to  harmless  waste  products  [1-4].  This  shuttling  of 
chemical  elements  between  their  oxidized  and  reduced  forms  is  known  as  nutrient  cycling  [5]. 
Nutrient  cycling  is  not  only  important  in  bioremediation,  but  is  also  an  essential  feature  of  the 
Earth’s  biogeochemical  cycles,  and  thus  forms  a  key  component  of  global  Earth-system  models 
[6-8].  Improving  our  approaches  to  bioremediation,  and  our  basic  understanding  of  how 
biogeochemical  cycles  work,  requires  us  to  be  able  to  predict  the  response  of  nutrient-cycling 
microbial  communities  to  environmental  changes  such  as  addition  of  pollutants.  Yet  despite  a 
wealth  of  new  information  emerging  from  DNA  sequencing  technologies  [9],  our  basic 
understanding  of  the  principles  underlying  the  structure  and  function  of  these  microbial 
communities  remains  very  poor  [10,1 1].  The  overarching  goal  of  our  project  was  to  reveal 
fundamental  principles  underlying  the  assembly  and  function  of  nutrient-cycling  microbial 
communities,  using  a  combination  of  mathematical  modeling  and  laboratory  experiments.  Our 
specific  aims  were: 

1.  To  characterize  experimentally  the  dynamical  changes  that  occur  during  the  development 
of  a  microbial  nutrient-cycling  ecosystem  and  to  develop  a  mathematical  model  that  can 
reproduce  these  changes. 

2.  Using  this  insight,  to  predict  how  these  systems  respond  to  environmental  change  and  to 
test  these  predictions  experimentally. 

3.  To  apply  these  principles  to  the  targeted  design  of  microbial  communities  for  applications 
in  bioremediation. 


To  achieve  these  aims,  our  projeet  eombined  theoretieal  and  experimental  work.  On  the 
theoretieal  side,  our  models  are  based  on  ordinary  differential  equations,  in  whieh  variables 
represent  the  population  densities  of  different  mierobial  metabolie  types,  and  the  eoneentrations 
of  key  ehemieal  speeies  [12,13].  The  growth  of  a  given  population  depends  on  the  coneentration 
of  its  food  sourees,  and  as  it  grows  it  exeretes  waste  produets  that  ean  be  used  as  food  by  other 
populations  [14,15].  On  the  experimental  side,  we  have  used  laboratory  mieroeosms,  made  from 
pond  sediment  and  water,  supplemented  by  earbon  and  sulfur  sourees  [16].  Over  several  months 
these  mieroeosms  develop  into  vertieally  stratified,  stable,  nutrient-eyeling  eommunities.  We  are 
able  to  traek  the  development  of  ehemieal  gradients  within  the  mieroeosms  using  voltammetry 
[17],  and  to  traek  the  mierobial  eommunity  eomposition  by  extraetion,  fingerprinting  and 
sequeneing  of  DNA  [18]. 

Summary  of  the  most  important  results 

Aim  1 :  Dynamieal  development  of  a  nutrient-eyeling  eeosystem 

Aim  1.1:  Experimental  characterization  of  microbial  ecosystem  development 

We  have  established  a  eomplete  set  of  ehemieal  data  (sulfide  and  oxygen  profiles)  and  mierobial 
eommunity  eomposition  data  (16S  metataxonomie  data  via  Illumina  sequeneing)  on  1 1 
mieroeosms  whieh  were  set  up  in  late  2013,  and  monitored  weekly  for  30  weeks  (Figure  1  and 
Figure  7).  The  mieroeosms  were  set  up  with  three  different  levels  of  added  eellulose,  whieh  we 
know  has  a  drastie  effeet  on  mieroeosm  development.  We  ineubated  these  mieroeosms  under 
fixed  light  and  temperature  eonditions  and  monitored  their  development  in  terms  of  ehemieal 
eomposition  (sulfide  and  oxygen  profiles)  on  a  weekly  basis  for  a  period  of  30  weeks.  We  also 
obtained  samples  for  DNA  analysis  of  the  mierobial  eommunity  eomposition,  non-destruetively 
using  a  thin  eorer.  As  the  mierobial  eommunity  develops,  sulfur-redueing  baeteria  eonvert  the 
sulfate  to  sulfide,  driven  by  aeetate  and  hydrogen  whieh  are  released  by  the  deeomposition  of  the 
eellulose.  Thus  traeking  the  produetion  of  sulfide  provides  a  measure  of  the  funetional  ehanges 
oeeurring  in  this  eeosystem.  Our  first  observation  was  that  the  total  amount  of  sulfide  in  our 
mieroeosms  peaks  at  a  time  and  level  whieh  depends  on  the  amount  of  earbon  souree  added.  We 
also  obtained  mueh  more  detailed  information  on  the  speeies  eomposition  of  our  mieroeosms,  by 
both  DGGE  fingerprinting  and  Illumina  sequeneing  of  samples  taken  non-destruetively  from  our 
mieroeosms.  The  results  show  that  speeies  eomposition  ehanges  strongly  in  time,  and  with  the 
initial  eellulose  and  sulfate  eoneentrations  (Figure  7).  They  also  reveal  large  blooms  of  Chlorobi 
(green  sulfur  baeteria)  in  the  mieroeosms  with  higher  initial  eellulose  eoneentrations,  eonfirming 
the  importanee  of  the  sulfur  eyele  in  these  eeosystem.  Plotting  the  sequenee  data  as  a  non-metrie 
multidimensional  sealing  (NMDS)  plot,  in  whieh  the  distanee  between  two  points  approximately 
represents  the  differenee  in  speeies  eomposition  between  two  samples  (Figure  7),  we  see  that 
mieroeosms  with  different  initial  eellulose  eoneentrations  follow  very  distinet  dynamieal 
trajeetories.  All  eommunities  ehange  initially  rapidly,  and  then  more  slowly,  over  the  15  week 
time  period.  The  greater  the  amount  of  initial  added  eellulose,  the  more  the  eommunity 
eomposition  ehanges  eompared  to  its  initial  inoeulating  sediment.  Combining  these  observations 
with  our  ehemieal  data,  whieh  shows  very  different  dynamieal  ehanges  in  sulfide  eoneentration 
depending  on  the  amount  of  initial  eellulose,  suggests  that  the  produetion  of  sulfide  by  sulfur 
redueing  mierobes,  driven  by  redueing  equivalents  generated  by  eellulose  degradation,  may 


ultimately  drive  community  development.  Intriguingly,  our  data  also  reveals  that  species 
diversity  correlates  inversely  with  sulfide  concentration,  initially  decreasing  as  our  microcosms 
develop,  but  later  increasing  again.  This  effect  is  most  pronounced  for  the  experiments  with  the 
largest  amount  of  added  cellulose;  suggesting  that  the  initial  diversity  loss  may  be  driven  by  the 
production  of  sulfide. 

We  have  also  performed  an  experiment  (revised  manuscript  submitted  to  Environmental 
Microbiology)  in  which  100  replicate  microcosms  were  set  up,  and  sampled  after  either  8  weeks’ 
or  16  weeks’  development.  These  microcosms  show  divergence  in  community  composition  as  a 
function  of  time  (which  we  do  not  see  in  larger  microcosms),  and  appear  to  form  alternative 
stable  states,  diverging  in  time  in  a  system  size-dependent  manner.  The  alternative  stable  states 
are  characterized  by  differences  in  the  degrader  communities;  one  group  of  steady  state 
microcosm  communities  is  enriched  with  specific  Firmicutes  species,  while  the  other  is  enriched 
with  specific  Bacteroidetes.  The  communities  dominated  by  members  of  these  two  phyla  also 
contain  distinct  populations  of  sulphate -reducing  bacteria.  This  biomodality  in  community 
composition  appears  to  arise  during  recovery  from  a  low-diversity  state  that  follows  initial 
cellulose  degradation  and  sulphate  reduction.  These  results  are  significant  because  relatively  few 
observations  of  alternative  stable  states  have  been  reported  for  microbial  ecosystems,  even 
though  they  are  rather  widely  reported  for  macro-ecosystems. 

Aim  1.2:  Mathematical  modeling  of  ecosystem  development 

During  the  course  of  this  project  we  have  developed  two  complementary  frameworks  for 
modelling  ecosystem  development  in  nutrient-cycling  microbial  communities.  Both  models  are 
based  on  ordinary  differential  equations  for  the  dynamics  of  the  various  microbial  populations 
and  chemical  species  in  our  ecosystems.  The  first  model  focuses  on  the  non-linearity  of 
microbial  growth  kinetics;  the  second  focuses  on  feedbacks  between  the  growth  of  different 
microbial  functional  types. 

Model  1:  non-linear  growth  (Appendix  1).  We  have  completed  an  extensive  theoretical  study  of 
perhaps  the  simplest  possible  mathematical  model  for  a  nutrient-cycling  ecosystem;  a  model  in 
which  two  populations  exchange  metabolites,  population  1  converting  the  reduced  form  of  a 
chemical  species  (which  we  denote  Sr)  into  its  oxidized  form  (so)  and  population  2  performing 
the  reverse  transformation.  This  model,  though  simple,  actually  provides  a  rather  generic 
representation  of  many  of  the  Earth’s  biogeochemical  cycles.  For  this  model,  the  differential 
equations  describing  the  growth  of  the  two  populations  and  the  concentrations  of  So  and  Sr  can  be 
solved  analytically  for  the  steady  state  (Appendix  1).  Importantly,  we  find  that  the  model 
undergoes  a  transition  between  a  “reduced”  ecosystem  state  (with  high  levels  of  Sr)  and  an 
“oxidized”  ecosystem  state  (with  high  levels  of  So)  as  we  vary  the  maximal  growth  rate  of  either 
the  reducing  or  the  oxidizing  microbial  population  (Figure  2).  We  termed  this  transition  a  “redox 
regime  shift”  in  analogy  with  sharp  transitions  in  ecosystem  state,  known  as  regime  shifts,  which 
have  been  observed  in  other  ecosystems  (eg  the  shift  between  temperate  and  desert  states  for 
savanna  ecosystems).  We  believe  that  this  transition  is  relevant  to  our  experiments,  because 
varying  the  maximal  growth  rate  of  the  reducer  population  provides  a  proxy  for  changes  in  the 
levels  of  acetate/hydrogen,  which  are  in  turn  produced  by  the  degradation  of  cellulose.  Thus, 
varying  the  concentration  of  cellulose  in  our  experiments  corresponds  to  varying  the  maximal 


growth  rate  parameter  in  our  very  simple  2-population  model.  The  advantage  of  using  sueh  a 
simple  model  is  that  we  ean  obtain  analytieal  solutions,  which  allow  us  to  access  the  fundamental 
mechanisms  underlying  the  regime  shifting  behavior.  We  find  that  the  regime  shift  occurs  when 
the  microbial  population  size  becomes  saturated,  so  that  a  population  can  no  longer  respond  to  a 
change  in  the  level  of  its  redox  substrate;  thus  these  regime  shifts  are  a  fundamental  consequence 
of  the  non-linearity  of  the  microbial  population  dynamics. 

To  verify  that  this  regime  shifting  behavior  is  actually  relevant  for  real  ecosystems,  we  also 
simulated  a  more  detailed,  2-box  model,  including  4  microbial  populations  (Figure  3).  In  this 
model,  cellulose  degraders,  which  release  acetate  from  cellulose,  and  photosynthesizers,  which 
produce  oxygen  when  stimulated  by  light,  are  explicitly  included.  This  model  mimics  spatial 
heterogeneity  in  a  basic  way  by  defining  separate  spatial  regions  representing  the  upper,  oxic 
zone  and  the  lower,  anoxic  zone  in  our  microcosms.  These  regions  are  coupled  by  chemical 
diffusion  and  within  each  region  we  assume  that  the  system  is  well-mixed.  We  have  shown  that 
this  model  also  predicts  a  drastic  change  in  ecosystem  state  in  response  to  gradual  changes  in 
environmental  parameters  (here  the  level  of  light  intensity  and  concentration  of  organic  matter). 
We  also  showed  that  the  parameters  for  which  such  redox  regime  shifts  are  predicted  by  the 
model  are  realistic  for  nutrient  cycling  populations  in  the  natural  environment. 

This  work,  which  has  been  published  in  Biogeosciences,  shows  that  proper  consideration  of  non- 
linearities  in  the  equations  modelling  microbial  population  dynamics  can  have  drastic  effects  on 
the  model  predictions.  This  is  an  important  conclusion  given  that  microbial  populations  in  large- 
scale  biogeochemical  models  are  usually  modelled  very  crudely. 

Model  2:  feedbacks  between  functional  groups  (Appendix  2)  We  have  also  collaborated  with 
Prof.  Jef  Huisman  and  Prof  Gerard  Muyzer  (University  of  Amsterdam,  The  Netherlands),  to 
develop  a  second  mathematical  model,  that  explains  how  alternative  stable  states  can  arise  in  a 
nutrient-cycling  microbial  ecosystem.  The  model  describes  the  coupled  population  dynamics  of 
cyanobacteria,  sulfate-reducing  bacteria  and  phototrophic  sulfur  bacteria,  as  found  in  our 
microcosms.  In  the  model,  cyanobacteria  produce  oxygen,  which  inhibits  the  sulfate-reducing 
bacteria  and  the  phototrophic  sulfur  bacteria.  Conversely,  the  sulfate -reducing  bacteria  produce 
sulfide,  which  inhibits  the  cyanobacteria.  The  phototrophic  bacteria  oxidize  sulfide,  creating  a 
sulfur  cycle.  Sulfur  can  also  be  oxidized  abiotically,  providing  an  alternative  cycling  mechanism 
(Figure  4).  Numerical  solution  of  this  model  shows  clearly  that,  in  some  regimes  of  parameter 
space,  the  model  has  two  alternate  stable  states,  one  of  which  is  oxic  and  rich  in  cyanobacteria, 
and  the  other  of  which  is  anoxic  and  rich  in  sulfate-reducing  bacteria.  The  anoxic  state  can  take 
two  forms,  depending  on  the  parameters;  a  form  in  which  phototrophic  bacteria  coexist  with 
sulfate-reducing  bacteria,  and  a  form  in  which  sulfide  is  oxidized  abiotically,  so  that  only  sulfate- 
reducing  bacteria  are  abundant.  Figure  5  shows  that  the  model  dynamics  leads  to  these  very 
different  final  states  when  we  make  a  small  change  in  the  initial  population  abundances. 

Aim  2:  Sensitivity  to  environmental  change 


Aim  2.1:  Predict  how  our  ecosystems  respond  to  changes  in  carbon,  sulfur  and  light 


Our  feedback-based  mathematical  model  sheds  new  light  on  how  nutrient-cycling  microbial 
ecosystems  may  respond  to  environmental  changes.  Specifically,  the  model  predicts  the 
existence  of  hysteresis:  changes  in  ecosystem  state  in  response  to  parameter  changes,  that  are  not 
reversible  when  the  parameter  change  is  reversed.  This  is  illustrated  in  Figure  6a,  where  we  show 
the  results  of  a  series  of  simulations  of  the  model,  where  in  each  successive  simulation  we 
increase  the  oxygen  diffusivity,  initiating  the  new  simulation  from  the  steady  state  of  the 
previous  simulation.  This  mimics  an  experiment  in  which  oxygen  diffusivity  is  slowly  increased 
(eg  by  changing  the  temperature).  Initially,  for  low  oxygen  diffusivity,  the  system  is  anoxic 
(oxygen  concentration  is  low).  As  oxygen  diffusivity  is  increased,  the  oxygen  concentration 
increases  slightly  (blue  line  in  Figure  6a),  before  making  a  sudden  transition  to  a  much  more  oxic 
state  (red  line  in  Figure  6a).  Starting  from  the  final,  oxic,  state  of  this  series  of  simulations,  we 
then  gradually  decrease  the  oxygen  diffusivity,  in  the  same  manner.  The  red  line  shows  the 
results:  oxygen  concentration  remains  high  over  a  range  of  oxygen  diffusivities  where  we 
obtained  the  anoxic  state  in  our  previous  set  of  simulations.  Eventually  we  do  see  a  transition 
back  to  the  anoxic  state. 

This  prediction  of  hysteresis  has  important  consequences  for  the  management  of  freshwater 
lakes.  For  example,  if  hysteresis  is  present  then  attempts  to  reoxygenate  anoxic  lakes  will  require 
much  more  oxygen  to  be  pumped  into  the  deep  water  than  a  simple  calculation  might  suggest. 
Interestingly,  we  have  evidence  that  hysteresis  does  occur  in  a  freshwater  lake  that  is  being 
studied  by  our  Dutch  collaborators:  Lake  Vechten  in  the  Netherlands.  This  lake  is  typically 
mixed  and  oxic  during  the  winter  (when  oxygen  diffusivity  is  high)  and  stratified  (with  an  upper 
oxic  layer)  during  the  summer,  when  oxygen  diffusivity  is  low.  Figure  6b  shows  a  plot  of  oxygen 
concentration  at  a  depth  of  7m  as  a  function  of  the  stratification  of  the  lake;  the  data  points  were 
sampled  at  different  times  over  the  course  of  a  year.  When  measurements  started  in  March  (blue 
data  point),  the  lake  was  mixed  and  oxic.  A  transition  to  a  stratified  state  occurred  over  the  next 
few  months,  and  the  oxygen  concentration  at  7m  dropped  dramatically.  However,  over  the  winter 
months,  when  the  lake  mixed  again,  it  remained  anoxic  for  several  weeks,  before  making  a 
sudden  transition  back  to  the  mixed  oxic  state.  Thus  the  lake  shows  clear  evidence  of  hysteresis. 
Metataxonomic  analysis  of  samples  from  the  lake  community  over  time  shows  negative 
correlations  in  the  abundance  of  cyanobacterial  with  sulfate-reducing  bacteria  and  phototrophic 
sulfur  bacteria,  while  the  two  latter  groups  are  positively  correlated;  also  in  agreement  with  the 
mathematical  model.  This  work  is  under  review  for  publication  in  Nature  Communications . 

Aim  2.2:  Compare  the  model  predictions  to  experimental  data 

As  explained  in  the  previous  section,  we  have  been  able,  at  least  qualitatively,  to  compare  the 
prediction  of  hysteresis  from  our  mathematical  model  to  field  data  from  a  freshwater  lake.  We 
have  also  been  carrying  out  microcosm  experiments  to  investigate  the  effects  of  environmental 
change  on  nutrient-cycling  microbial  communities,  as  detailed  in  our  previous  reports. 

We  have  also  performed  microcosm  experiments  in  which  we  have  obtained  and  analyzed  a 
complete  set  of  data  for  a  series  of  microcosms  containing  different  concentrations  of  cellulose 
and  sulfate.  After  incubating  these  microcosms  for  16  weeks,  we  measured  vertical  profiles  of 
sulfide  and  oxygen,  and  took  samples  for  DNA  analysis  (via  Illumina  16S  rRNA  gene 
sequencing).  We  also  analyzed  the  DNA  sequence  data  in  terms  of  microbial  functional  groups. 


finding  that  different  funetional  groups  showed  different  responses  to  environmental  ehange:  for 
example  sulfate  reducers  maintain  stable  abundance  and  diversity  in  response  to  changes  in 
sulfate  and  cellulose,  but  oxygenic  phototrophs  show  a  collapse  in  diversity  as  cellulose 
concentration  increases.  We  believe  that  these  results  may  be  useful  in  understanding  the  niche 
structure  of  microbial  communities.  Within  a  particular  microbial  functional  group,  we  suppose 
that  individual  taxa  occupy  specific  niches  (for  example  defined  by  a  particular  pH  range,  spatial 
location  etc).  In  our  experiment,  we  can  track  how  the  occupation  of  these  niches  changes  as  we 
change  the  environmental  conditions.  For  example  for  sulfate-reducing  bacteria,  we  see  that 
diversity  is  higher  in  microcosms  where  the  sulfate-reducing  population  is  higher  (diversity  is 
positively  correlated  with  abundance);  Figure  8a.  This  suggests  that  more  niches  are  present  in 
microcosms  with  higher  sulfate  reducer  diversity.  However,  we  observe  the  opposite  for 
phototrophic  sulfur  bacteria:  diversity  is  actually  lower  in  microcosms  where  phototrophic  sulfur 
bacteria  are  more  abundant;  Figure  8b.  This  can  be  explained  by  a  picture  in  which  a  few  niches 
expand  in  size  as  the  total  population  expands,  but  most  niches  remain  the  same.  For  oxygenic 
phototrophs,  we  see  yet  another  scenario:  diversity  remains  the  same  as  the  abundance  increases; 
Figure  8c.  This  is  consistent  with  a  picture  in  which  the  number  of  niches  remains  the  same  but 
they  all  increase  in  size  proportionally.  To  our  knowledge,  this  is  the  first  observation  to  suggest 
that  the  ecology  of  different  microbial  functional  groups  might  be  fundamentally  different  - 
suggesting  that  attempts  to  apply  “ecological  theory”  as  defined  for  macro-organisms,  to 
microbial  communities,  may  well  be  over-simplistic.  We  are  working  on  a  manuscript  describing 
these  results. 

Aim  3:  Targeted  optimization  of  microbial  communities  for  bioremediation 

We  have  established  a  model  system  for  studying  the  biodegradation  of  the  pollutant  PCP  in  our 
microcosms,  using  the  biodegrading  organism  D.  hafniense,  and  assays  for  both  D.  hafniense  (by 
PCR)  and  PCP  (using  gas  chromatography-mass  spectrometry  (GC-MS)).  This  system  is 
interesting  because  D.  hafniense  requires  hydrogen  as  an  electron  donor.  The  sulfate  reducers, 
which  play  a  key  role  in  our  microcosms,  may  inhibit  D.  hafniense,  by  competing  for  hydrogen 
(which  they  use  to  reduce  sulfate).  However,  some  literature  suggests  that  sulfate  reducers  may 
under  some  circumstances  actually  aid  the  growth  of  D.  hafniense  by  producing  hydrogen  (which 
happens  when  they  grow  fermentatively).  Motivated  by  this  potentially  interesting  conflict,  we 
have  set  up  a  range  of  microcosms  with  different  amounts  of  added  sulfate,  some  of  which  have 
been  supplemented  by  PCP  and  /  or  D.  hafniense.  From  a  practical  point  of  view,  an  especially 
interesting  question  concerns  the  length  of  time  that  a  population  of  D.  hafniense  is  maintained  in 
the  community.  Previous  studies  of  bioaugmentation  (in  which  a  pollutant-degrading  organism  is 
added  to  a  microbial  community)  suggest  that  the  added  population  tends  to  die  out  within  a  few 
days,  limiting  the  usefulness  of  the  method  for  pollutant  degradation.  Our  data,  in  contrast, 
suggest  that  D.  hafniense  can  persist  in  our  microcosms  for  many  weeks. 

Aim  3.1:  Theoretical  predictions  of  optimal  conditions  for  PCP  degradation 
We  have  extended  our  “feedback”  mathematical  model  (Figure  4),  to  include  a  population  of 
reductive  dehalogenators  (representing  D.  hafniense).  This  population  consumes  PCP,  competes 
with  sulfate-reducing  bacteria  for  hydrogen,  and  competes  with  all  other  organisms  for 
phosphorous.  It  is  also  inhibited  by  oxygen  (since  D.  hafniense  is  anaerobic).  The  extended 
model  structure  is  shown  in  Figure  9.  Our  model  predicts  a  strong  dependence  of  the  ecosystem 


state  on  the  total  amount  of  sulfur  that  is  present.  At  low  total  sulfur,  cyanobacteria  dominate  and 
the  system  is  oxic,  with  no  PCP  degradation.  At  higher  levels  of  sulfur,  the  action  of  sulfate- 
reducing  bacteria  causes  the  system  to  become  anoxic.  If  the  sulfur  level  is  not  too  high  then 
sulfate  reducers  are  eventually  outcompeted  by  reductive  dehalogenators,  due  to  the  competition 
for  hydrogen.  However  if  the  level  of  sulfur  is  very  high  then  the  sulfate  reducers  will 
outcompete  the  reductive  dehalogenators,  leading  to  a  final  state  with  a  sulfur  cycle  but  without 
PCP  degradation.  Thus,  the  model  predicts  that  PCP  degradation,  and  persistence  of  D. 
hafniense,  will  only  occur  at  intermediate  concentrations  of  added  sulfur. 

Aim  3.2:  Experimental  measurements  of  PCP  degradation 

We  have  carried  out  experiments  using  microcosms  supplemented  with  PCP  and  D.  hafniense, 
compared  to  unsupplemented  microcosms,  in  the  presence  of  different  levels  of  added  sulfate. 
Our  results  show  the  supplementation  with  PCP  affects  community  structure  drastically,  and, 
excitingly,  that  D.  hafniense  can  persist  in  our  microcosms  over  long  times  (Figure  10).  This  is  in 
contrast  to  many  attempts  at  bioaugmentation,  including  with  D.  hafniense  in  soil  microcosms,  in 
which  the  added  organism  is  rapidly  lost.  Consistent  with  our  model,  we  see  that  D.  hafniense  is 
lost  more  rapidly  from  microcosms  without  added  sulfate  than  from  microcosms  to  which  a  low 
concentration  of  sulfate  has  been  added.  We  also  see  strong  changes  in  microbial  community 
composition  depending  on  the  levels  of  added  PCP  and  sulfate.  We  have  also  just  obtained 
sequence  data  from  a  final  experiment,  completed  last  year,  in  which  a  higher  level  of  sulfate 
was  added  to  PCP-  and  D.  /zq/h/en^e-supplemented  microcosms.  With  the  assistance  of  a  Masters 
student,  we  will  be  analyzing  these  data  over  the  next  few  months  to  determine  whether,  as 
predicted,  the  high  level  of  sulfate  allows  the  sulfate  reducers  to  out-compete  D.  hafniense  and 
shorten  its  persistence  in  the  microcosm  system. 

Taken  together,  our  modelling  and  simulation  results  on  the  PCP-O.  hafniense  microcosm 
system  provide  new  insights  into  the  mechanisms  that  may  come  into  play  when  reductive 
dehalogenators  are  added  to  a  nutrient-cycling  microbial  ecosystem.  It  is  clear  from  our  results 
that  understanding  the  structure  and  function  of  nutrient-cycling  microbial  communities  can 
indeed  lead  to  targeted  strategies  for  bioremediation;  producing  much  longer  persistence  of 
added  pollutant-degrading  organisms  than  has  been  achieved  before. 
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Figure  1 

Photograph  of  the  microcosms  which  we  have  been  tracking  dynamically  as  they  undergo  community  assembly. 
This  photograph  was  taken  10  weeks  into  the  experiment.  The  different  colored  labels  indicate  microcosms  with 
different  added  amounts  of  cellulose. 
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Figure  2 

Redox  regime  shifts  in  model  nutrient  cycles,  obtained  by  analytical  solution  of  a  simple  model.  In  this  model,  in 
the  “fully  biotic  case” ,  one  population  transforms  a  chemical  from  its  oxic  form  So  to  its  anoxic  form  and  another 
performs  the  reverse  transformation.  In  the  “biotic-abiotic  case” ,  a  microbial  population  performs  the  oxic  to  anoxic 
transformation  but  the  anoxic  to  oxic  transformation  is  mediated  by  an  abiotic  reaction.  The  global  redox  state,  as 
measured  by  the  oxidized  fraction  So/stot,  predicted  by  the  steady-state  solution  of  the  model  equations  for  the  fully 
biotic  cycle  (a  and  c)  or  the  biotic-abiotic  cycle  (b  and  d)  is  plotted  as  a  function  of  parameters  that  form  proxies  for 
the  degree  of  reductive  or  oxidative  driving.  These  parameters  are:  for  reductive  driving,  the  maximal  growth  rate  of 
the  reductive  population,  Vot  (a  and  b),  and,  for  oxidative  driving,  either  the  maximal  growth  rate  of  the  oxidative 
population  v^o  (c)  or  the  maximal  abiotic  oxidation  rate  Ua  (d).  The  results  show  a  shift  between  oxidized  and  reduced 
ecosystem  states  as  a  threshold  in  reductive  or  oxidative  driving  is  crossed;  the  sharpness  of  this  transition  increases 
with  the  concentration  of  the  chemical  species  being  cycled,  Stot  (in  panels  a  to  c,  red  to  blue  lines;  red:  Stot  =  20/rM, 
pink:  stot  =  0.2mM,  purple:  Stot  =  2mM  blue:  Stot  =  20mM;  in  panel  d,  black  line  Stot  =  2M:  dark  purple  line: 
Stot  =  200mM). 
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Figure  3 

The  “four-population,  two-box”  mathematical  model  and  its  predictions  for  redox  regime  shifts,  (a)  Illustration 
of  the  model.  Oxidative  and  reductive  processes  take  place  in  separate  spatial  zones,  linked  by  chemical  diffusion. 
The  model  explicitly  represents  the  population  dynamics  of  microbial  photosynthesizers,  decomposers,  reducers  and 
oxidizers,  and  the  chemical  dynamics  of  oxygen,  So,  Sr  and  acetate.  Light  intensity  and  organic  matter  availability 
are  treated  as  control  parameters.  The  dynamical  equations  corresponding  to  the  model  are  integrated  numerically 
to  find  the  steady-state  solution,  (b)  Steady-state  solution  of  the  complete  ecosystem  model,  obtained  numberically, 
plotted  as  a  function  of  the  control  parameters,  light  intensity  (relative  to  the  typical  value  10  /iEinstein  s“^m“^) 
and  organic  matter  concentration  (relative  to  the  typical  value  100  mg  cm“^).  The  color  represents  the  global  redox 
state  (see  color  key).  The  model  shows  redox  regime  shifts  as  the  organic  matter  concentration  is  varied  at  fixed  light 
intensity  (vertical  dashed  line)  or  as  the  light  intensity  is  varied  at  fixed  organic  matter  concentration  (horizontal 
dashed  line). 


A  schematic  diagram  of  our  microbial  ecosystem  model  which  shows  alternate  stable  states.  Diamond  shaped  boxes 
represent  microbial  functional  groups  and  circles  represent  substrates.  CB:  Cyanobacteria,  PB:  Phototrophic  sulfur 
bacteria,  SB:  Sulfate  reducing  bacteria,  P:  Phosphorus,  O:  Oxygen,  S^:  Reduced  sulfur  and  Sq-  Oxidized  sulfur. 
Dark  blue  arrows  denote  consumption  of  a  chemical  by  the  microbial  population  that  they  point  to.  Red  arrows 
denote  production  of  a  chemical  by  the  microbial  population  that  they  originate  from.  Orange  lines  from  a  chemical 
to  a  microbial  population  represent  growth  inhibition  of  the  microbial  population  by  that  substrate.  The  green  arrow 
denotes  the  abiotic  reaction  of  reduced  sulfur  and  oxygen  to  produce  oxidized  sulfur. 
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Evidence  for  alternate  stable  states  in  our  mathematical  model  represented  by  sensitivity  to  the  initial  cyanobacterial 
population  density.  We  use  a  log  scale  for  the  time  axis  to  visualize  the  fact  that  different  processes  in  the  system 
occur  on  different  timescales.  For  all  plots  initial  npB  =  nsB  =  1  x  10^^  and  So,b  =  Sp^B  =  500/rM,  Pb  =  O.bfiM, 
Ob  =  300.0/rM.  (a)  shows  the  dynamics  of  the  model  for  an  initial  ncB  =  1  x  10®.  (b)  shows  the  dynamics  of 
the  model  for  an  initial  ncB  =  1  x  10®.  The  plot  shows  population  densities  as  a  function  of  time:  Green  =  ncB  j 
Orange  =  npB  and  Purple  =  nsB  ■  In  (a)  the  model  predicts  a  final  state  with  coexisting  sulfate-reducing  bacteria 
dn  phototrophic  sulfur  bacteria;  in  (b)  it  predicts  a  final  state  with  only  cyanobacteria. 


Figure  6 

Hysteresis  in  a  nutrient-cycling  microbial  community,  (a)  Predictions  of  our  mathematical  model.  The  results  show 
the  effect  of  increasing  and  then  decreasing  the  oxygen  diffusivity  (see  orange  arrows)  where  the  initial  state  of  each 
simulation  is  set  to  the  final  microbial  ecosystem  state  of  the  previous  simulation.  The  forward  and  reverse  trajectories 
do  not  follow  the  same  path,  (b)  Evidence  for  regime  shifts  between  oxic  and  anoxic  states  in  lake  Vechten.  Hysteresis 
in  the  strength  of  lake  stratification  as  a  function  of  oxygen  saturation  at  7m  depth,  from  March  2013  (blue  square) 
to  April  2014  (red  square),  following  the  red  arrows.  Again  the  forward  and  reverse  trajectories  follow  different  paths; 
the  system  becomes  temporarily  stuck  in  a  mixed  anoxic  state. 
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Figure  7 

Community  composition  of  our  microcosms  as  a  function  of  time,  during  the  process  of  community  assembly,  as 
assessed  by  16S  rRNA  gene  sequencing.  Microcosms  have  been  sampled  weekly,  for  3  different  initial  concentrations 
of  added  cellulose,  here  denoted  100%,  45%  and  15%  (relative  to  a  standard),  (a):  Community  composition  at  the 
phylum  level.  Different  colors  represent  different  phyla.  The  left  2  bars  show  the  initial  pond  sediment  used  to 
create  the  microcosms,  (b)  The  same  data,  represented  on  an  NMDS  plot,  where  the  distance  between  two  points 
approximates  the  difference  in  composition  between  two  samples.  It  can  clearly  be  seen  that  microcosms  with  different 
initial  amounts  of  cellulose  follow  different  trajectories  of  community  composition,  (c)  Community  species  diversity, 
as  measured  by  the  number  of  observed  taxa,  as  a  function  of  time,  for  three  microcosms  with  differing  amounts  of 
added  cellulose. 
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Figure  8 

Results  from  our  microcosm  experiments  with  different  amounts  of  added  cellulose  and  sulfate.  For  three  functional 
groups  of  bacteria  (as  defined  by  taxonomic  classification),  we  plot  the  diversity  (as  measured  by  Shannon  Index)  as 
a  function  of  the  relative  abundance  of  individual  taxa  (determined  at  the  97%  level) .  The  three  microbial  functional 
groups  show  very  different  trends  in  diversity-abundance  relation. 


Schematic  illustration  of  our  mathematical  model  for  a  microbial  nutrient-cycling  ecosystem  supplemented  by 
reductive  dehalogenators.  Symbols  are  as  in  Figure  4.  Haf  represents  the  population  of  D.  Hafniense,  H  represents 
hydrogen  and  PCP  represents  the  pollutant  PCP. 
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Figure  10 

Results  of  our  bioremediation  experiment.  In  this  experiment,  microcosms  were  set  up  with  and  without  addition 
of  the  chemical  pollutant  PCP,  with  and  without  addition  of  sulfate.  Microcosms  were  also  supplemented  with  the 
PCP-degrading  organism  B.  hafniense  PCP-1,  either  at  the  start  only  (DHl  in  the  plot)  or  both  at  the  start  and 
after  1  month  (DH2  in  the  plot).  The  abundance  of  D.  hafniense  PCP-1  was  monitored  by  RTPCR  at  different  times, 
with  the  results  shown  by  the  colored  bars.  Our  results  show  that  in  microcosms  that  are  not  supplemented  by  sulfate 
(left  3  datasets),  D.  hafniense  PCP-1  dies  out  over  time.  However,  in  microcosms  that  are  supplemented  by  sulfate 
(right  2  datasets),  D.  hafniense  PCP-1  is  maintained  in  the  community  over  much  longer  times.  This  suggests  that 
understanding  nutrient  cycling  can  provide  strategies  for  successful  bioremediation  via  bioaugmentation. 


Appendix  1:  The  four-population,  two-box  model 


Here  we  list  the  equations  corresponding  to  the  four-population,  two-box  model  illustrated  in  Fig.  3. 


dnp{t) 

dt 


vpnp{t)L 
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Eq.  (1)  describes  growth  of  the  photosynthesizer  population,  with  density  np.  Photosynthesizers  are  located  in 
the  upper  box.  They  grow  in  response  to  light  intensity  L  (which  is  a  constant  parameter  in  our  model),  with 
maximal  growth  rate  vp  and  half-saturation  constant  Kp.  We  impose  a  population  size  limitation  term  with  maximal 
population  size  np^max  and  a  death  rate  d. 
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Eq.  (2)  describes  growth  of  the  decomposer  population,  with  density  np,.  Decomposers  are  located  in  the  lower 
box.  They  grow  in  response  to  organic  carbon  concentration  C  (which  is  a  constant  parameter  in  our  model),  with 
maximal  growth  rate  vp,  and  half-saturation  constant  Kq.  We  impose  a  population  size  limitation  term  with  maximal 
population  size  nD,max  and  a  death  rate  d. 
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Eq.  (3)  describes  growth  of  the  microbial  reducer  population,  with  density  Uor-  Reducers  are  located  in  the  lower 
box.  Reducer  growth  requires  both  the  oxidized  form  of  the  redox  chemical  species,  So,  and  acetate.  We  describe  this 
using  multiplicative  Monod  terms,  with  maximal  growth  rate  Wor  and  half-saturation  constants  for  Sq  and 
for  acetate.  The  concentration  of  acetate  is  denoted  as  a{t).  The  reducer  population  size  is  limited  to  nor, max  and  we 
impose  a  death  rate  d. 
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Eq.  (4)  describes  growth  of  the  microbial  oxidizer  population,  with  density  n^o-  Oxidizers  are  located  in  the  upper 
box.  Oxidizer  growth  requires  both  the  reduced  form  of  the  redox  chemical  species,  Sr,  and  oxygen.  Again,  we  describe 
this  using  multiplicative  Monod  terms,  with  maximal  growth  rate  Vro  and  half-saturation  constants  ATro  for  Sr  and 
ATox  for  oxygen.  The  concentration  of  oxygen  is  denoted  as  o{t).  The  oxidizer  population  size  is  limited  to  nro.max 
and  we  impose  a  death  rate  d. 
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Eq.  (5)  describes  the  dynamics  of  the  oxygen  concentration.  Oxygen  is  assumed  to  be  located  only  in  the  upper 
box.  Oxygen  is  produced  upon  growth  of  the  photosynthesizers  (first  term  in  Eq.  (5)),  with  yield  parameter  7ox 
(number  of  micromoles  of  oxygen  produced  per  bacterial  division  cycle).  Oxygen  is  consumed  upon  growth  of  the 
oxidizer  population  (second  term  in  Eq.  (5)),  with  yield  parameter  7.  We  also  include  a  term  describing  consumption 
of  oxygen  by  other  processes,  not  included  in  the  model  (third  term  in  Eq.  (5)),  such  as  growth  of  aerobes  or  abiotic 
oxidation  reactions.  The  parameter  /3ox  controls  the  strength  of  this  “competition”  term. 
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Eq.  (6)  describes  the  dynamics  of  the  acetate  concentration.  Acetate  is  assumed  to  be  located  only  in  the  lower  box. 
Acetate  is  produced  upon  growth  of  the  decomposers  (first  term  in  Eq.  (6)),  with  yield  parameter  7ac.  Acetate  is 
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consumed  upon  growth  of  the  reducer  population  (second  term  in  Eq.  (6)),  with  yield  parameter  7  (here  assumed 
to  be  the  same  as  that  of  the  reducers).  We  also  include  a  “competition”  term  describing  consumption  of  acetate  by 
other  processes  (third  term  in  Eq.  (6)),  such  as  growth  of  methanogens.  The  parameter  /3ac  controls  the  strength  of 
this  term. 
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Eqs.  (7)-(10)  describe  the  dynamics  of  the  reduced  and  oxidized  forms  of  the  redox  chemical  species,  Sr  and  So,  in 
the  upper  and  lower  boxes. 


These  equations  were  solved  numerically  for  the  steady  state,  using  Euler  integration,  to  obtain  the  data  shown  in 
Fig.  3. 


Appendix  2:  Description  of  the  mathematical  model  for  a  nutrient-cycling  ecosystem  showing  alternate 

stable  states 

We  model  a  microbial  ecosystem  containing  three  microbial  functional  groups:  oxygen-producing  cyanobacteria 
(CB),  phototrophic  sulfur  bacteria  (PB)  such  as  purple  or  green  sulfur  bacteria  that  oxidize  reduced  sulfur  and 
produce  oxidized  sulfur,  and  sulfate-reducing  bacteria  (SB)  that  reduce  oxidized  sulfur  and  produce  reduced  sulfur. 
The  dynamics  of  our  mathematical  model  are  controlled  by  feedbacks  between  the  three  functional  groups.  The  oxygen 
produced  by  the  cyanobacteria  is  assumed  to  be  inhibitory  to  both  the  sulfate-reducing  bacteria  and  phototrophic 
sulfur  bacteria.  Conversely,  the  reduced  sulfur  produced  by  the  sulfate-reducing  bacteria  is  assumed  to  be  toxic  to 
cyanobacteria.  Reduced  sulfur  is  oxidized  by  the  phototrophic  sulfur  bacteria,  and  oxidized  sulfur  is  reduced  by 
the  sulfate-reducing  bacteria,  connecting  these  species  in  a  nutrient  cycle  as  they  pass  the  sulfur  back  and  forth. 
Importantly,  the  oxidation  of  reduced  sulfur  can  also  occur  abiotically,  producing  oxidized  sulfur  in  the  process. 

The  dynamics  of  this  model  are  described  by  Eqs.  1-7.  Eqs.  1-3  describe  the  dynamics  of  the  cyanobacteria, 
phototrophic  sulfur  bacteria  and  sulfate- redncing  bacteria  respectively,  with  population  densities  represented  by  ncB  j 
npB  and  nsB  respectively.  For  all  functional  groups  a  constant  removal  rate  D  is  applied  (e.g.  due  to  grazing  or 
viral  lysis).  Gi  and  It  represent  growth  and  inhibition  functions  respectively  for  population  i.  The  growth  of  each 
population  is  assumed  to  depend  on  the  availability  of  phosphorus,  a  key  limiting  resource  of  microbial  ecosystems  in 
many  aquatic  environments.  Additionally,  the  growth  of  the  phototrophic  sulfur  bacteria  depends  on  the  availability 
of  reduced  sulfur  and  the  growth  of  sulfate-reducing  bacteria  depends  on  the  availability  of  oxidized  sulfur. 


=  ^^CBGcB(R)7cB(-S'R)  7?ncB 
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Eqs.  4-7  describe  the  dynamics  of  the  substrates:  oxidized  sulfur,  reduced  sulfur,  oxygen  and  phosphorus  respec¬ 
tively.  So  represents  the  concentration  of  oxidized  sulfur  and  5'r  represents  the  concentration  of  reduced  sulfur.  O 
represents  the  concentration  of  oxygen  and  P  represents  the  concentration  of  freely  available  phosphorus.  For  every 
chemical  substrate  a  background  concentration  is  imposed  (labelled  by  subscript  b).  The  substrate  concentration 
tends  to  be  restored  to  this  background  concentration  at  a  rate  that  depends  on  the  difference  between  the  current 
concentration  and  the  background  concentration,  and  on  a  parameter  a.  The  parameter  a  is  specific  to  each  substrate, 
for  example  ao  controls  the  rate  at  which  oxygen  enters  or  leaves  the  system  (which  in  the  environment  could  repre¬ 
sent  seasonal  temperature-related  changes  in  oxygen  diffusivity) .  Importantly,  we  also  include  the  chemical  oxidation 
of  reduced  sulfur,  which  proceeds  at  a  rate  R  =  — A:2<S'rO,  consuming  oxygen  (O)  and  producing  oxidized  sulfur  (S'o) 
in  the  process. 
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One  could  in  principle  choose  many  different  forms  for  the  growth  and  inhibition  functions  Gi  and  1^.  We  expect  our 
qualitative  results  to  be  largely  insensitive  to  the  detailed  choice  of  growth  and  inhibition  function.  In  this  work,  we 
have  chosen  to  use  Monod  functions.  Eqs.  8-10  give  the  growth  functions  for  the  cyanobacteria  (Gcb),  phototrophic 
sulfur  bacteria  (Gpb)  and  sulfate  reducing  bacteria  (Gsb)  respectively. 
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In  Eqs.  8-10,  the  microbial  growth  rate  as  a  function  of  substrate  concentration  is  described  by  a  Monod  function 
of  the  form  where  v  (hr“^)  is  the  maximal  growth  rate  and  Kj  (/iM)  is  the  concentration  of  substrate  j  at 

which  the  growth  rate  is  half  the  maximal  growth  rate.  For  functional  groups  that  are  limited  by  the  availability 
of  more  than  one  nutrient  (the  phototrophic  sulfur  bacteria  and  the  sulfate-reducing  bacteria)  multiplicative  Monod 

kinetics  of  the  form  v  ^  used. 

Eqs.  11-13  describe  the  inhibition  functions  for  the  cyanobacteria  (Icb),  phototrophic  sulfur  bacteria  (/pb)  and 
sulfate- reducing  bacteria  (/sb)  respectively.  The  inhibition  of  microbial  growth  is  described  by  a  version  of  the 
well-known  Haldane  equation  ^  ^  . 
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